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Abstract 



A phenomenological model of the evolution of an ensemble of interact- 
ing dislocations in an isotropic elastic medium is formulated. The line-defect 
microstructure is described in terms of a spatially coarse-grained order param- 
eter, the dislocation density tensor. The tensor field satisfies a conservation 
law that derives from the conservation of Burgers vector. Dislocation motion 
is entirely dissipative and is assumed to be locally driven by the minimization 
of plastic free energy. We first outline the method and resulting equations of 
motion to linear order in the dislocation density tensor, obtain various sta- 
tionary solutions, and give their geometric interpretation. The coupling of 
the dislocation density to an externally imposed stress field is also addressed, 
as well as the impact of the field on the stationary solutions. 
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I. INTRODUCTION 



The presence and motion of dislocations in materials is known to affect many of their 
mechanical properties, including, for example, resistance to deformation and plastic response 
to relatively large stresses [|l|. Given the impact of dislocation dynamics on material proper- 
ties, there have been a number of models advanced attempting to describe the motion of line 
defects in response to both internal and applied stresses, which have substantially furthered 
our understanding of these defect-property correlations 

The starting point of a large class of models of dislocation dynamics is the calculation of 
direct interactions between a dislocation and a stress field, as given by the Peach-Koehler 
force P|. On the one hand it is possible to construct a Newtonian model of dislocation 
motion by combining the law of conservation of Burgers vector with Newton's second law, 
a Hookean stress-strain constitutive law and the differential relation between plastic strain 
and dislocation density, as has been done elsewhere. While this approach has been very 
useful in elucidating some important features of dislocation motion, its inherent limitations 
to relatively short time scales and small numbers of dislocations make its generalization to 
systems with large numbers of dislocations which evolve over long periods of time somewhat 
impractical. Numerical calculations of the dynamical evolution of an ensemble of pairwise 
interacting dislocations have been carried out with Molecular Dynamics methods to study 
relatively short time scales of the order of several vibrational periods A different approach 
to dislocation dynamics is based on purely dissipative or Langevin dynamics |^. While the 
level of description is still microscopic, the motion of each dislocation line is assumed to 
be overdamped. This method allows the study of longer time scales at the price of added 
phenomenology in the introduction of friction or damping forces. 

While the aforementioned approaches have proved quite fruitful, they share several draw- 
backs that limit their usefulness in studying the evolution of a large ensemble of disloca- 
tions over relatively long periods of time. Among them, we note the long ranged nature of 
dislocation-dislocation interactions which must be truncated in actual calculations, often in 
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some ad- hoc fashion In addition, the microscopic nature of the description places a se- 
vere restriction to the number of dislocation lines that can be considered. As a consequence, 
the scale of the resulting microstructure is essentially set by the number of dislocations 
in a typical simulation cell. We note, however, that recent computational advances have 
permitted studies involving up to 10® dislocations [0 . 

Experimental evidence, on the other hand, suggests that, at least in some cases, the 
mechanical response of a material may depend only weakly on small-scale microstructural 
details, and more strongly on the macroscopic density of dislocations. Such evidence 
implies that an approach based, a priori, on the dislocation density might best capture 
the important macroscopic features of such systems. In this paper we follow this approach 
and formulate a simplified, linear phenomenological model of dislocation motion, in terms 
of the dislocation density tensor. The model properly incorporates necessary conservation 
laws related to the topological structure of the dislocation lines, as well as the true tensorial 
character of the density. By contrast with previous approaches to this problem, and in the 
same spirit as irreversible thermodynamics, we employ a non-deterministic description which 
focuses on a small number of macro-variables and treats the remaining degrees of freedom, 
which are presumed to evolve on a faster time scale, as an effective heat bath. This is formally 
accomplished by applying a coarse-graining procedure to a system of dislocations in an 
elastically isotropic medium, and selecting the dislocation density tensor as the appropriate 
order parameter to describe the system at this scale. Of course, this approach is formally 
equivalent to Langevin type models mentioned earlier: dislocation motion is overdamped 
because of its interaction with other microscopic degrees of freedom that act as a "heat 
bath" . We argue, however, that the overdamped character of dislocation motion arises only 
at a semi-macroscopic scale, and that it includes not only interactions with the "heat bath" , 
but also mutual dislocation interactions within a coarse-graining cell. 

Our approach also parallels the study of the equilibrium statistical mechanics of a system 
of unbound dislocations undertaken by Nelson and Toner They investigated the effect 
of a collection of dislocations on the broken translational and rotational symmetries of 
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the crystalline state by examining the Fourier transform of the dislocation density- density 
correlation function at long wavelengths. In so doing, they argued that dislocations destroy 
long-ranged translational order, but not long-ranged orientational order. Central to this 
calculation is the identification of the appropriate Boltzmann weight (and therefore the free 
energy) that determines the probability of the occurrence of a given fluctuation, and of the 
corresponding defect configuration. In their view the relevant dislocation driving forces can 
be constructed from a free energy which is a functional of the dislocation density tensor, with 
components p^j, and a non-singular strain, with components efj^ , which describes vibrational 
degrees of freedom in the crystal. The contribution to the free energy that depends on pij 
was chosen equal to the energy of interaction between dislocations, as derived by Kosevich 
T0| . The logical next step, which we propose here, is to use the minimization of the singular 



contribution to the free energy as the driving force for dislocation motion. 

We finally mention complementary earlier work by Brand and Kawasaki [p!T| -[T3[], who 
derived phase equations for generic systems that support topological defects. A phase equa- 
tion involves a phase variable that changes slowly on the scale of the lattice spacing and 
in time, and that has a nonzero winding number around a topological defect. This type 
of equation represents the long-wavelength distortion of some ideal basic pattern and has 
been successfully used in a large variety of nonequilibrium systems fl^. We focus our at- 
tention on the case of a crystalline solid and retain the full dislocation density tensor in our 
description, since we believe that this is the natural coarse-grained variable to describe a de- 
fected crystalline solid. Our choice of variables may prove useful in two different directions. 
First, the equations that we present below may provide a starting point for direct numerical 
calculations of dislocation motion, a subject of considerable interest. Second, coupling to 
other slow variables of interest can be accomplished quite naturally. These would include, 
for example, energy or mass densities. The latter are necessary to describe other irreversible 
phenomena such as the effect of mass diffusion on dislocation motion. 

This paper is organized as follows. In Sec. ^ we describe the dynamical equations 
governing the evolution of the dislocation density in an isotropic elastic medium. In Sec. 



PI stationary solutions are identified and given geometrical interpretations, both in systems 
with and without externally imposed stresses. Finally, Sec. |IV| contains a short summary 
and discussion of future work. 

II. MODEL EQUATIONS 

Consider an isothermal, three-dimensional, linear isotropic elastic medium which contains 
a large, finite number of unbound dislocation lines which are free to move throughout the 
system. Suppose further that this medium is subdivided into a number of hypothetical 
averaging cells, each of volume Vq, such that a relatively large number of dislocation lines 
thread a given volume element. Now, if one wishes to describe the defect properties of this 
system on a length scale larger than (VoY^^, then one can define a coarse-grained variable, the 
dislocation density, which reflects the average dislocation content of a cell. As a dislocation is 
characterized by two directions in space, its line direction, x, and the direction of its Burgers 
vector, b, it is expected that the dislocation density must also carry such information and 
would, therefore, be a second-rank tensor. The components of this tensor will be denoted 
by pij hereafter. 

In modeling the dynamics of the defected system described above at this coarse-grained 
level, it is necessary to include "slow" variables, the relaxation of which takes place on time 
scales much longer than microscopic times p!5|JT6|| . Variables that satisfy global conserva- 
tion laws or that arise as a consequence of a spontaneously broken symmetry satisfy this 
requirement. Although line defects are not present in equilibrium in a three dimensional 
crystalline solid, we take the view that there exists a time scale, slow compared to rele- 
vant experimental times, in which the net Burgers vector in a crystalline solid is conserved. 
In practice, of course, dislocations may migrate to the surface of the crystal or to a grain 
boundary and disappear. Conservation of Burgers vector can be expressed in terms of the 
dislocation density tensor as pO[] 
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where eum are the components of the Levi-Civita tensor density, jmj are the components of 
a dislocation flux density tensor which govern the flow of Burgers vector in the medium and 
we have adopted the convention that repeated indices are summed over. The dislocation 
density tensor satisfies the relation, 



where the integral extends over an arbitrary surface in the crystal, and bj is the j-component 
of the sum of Burgers vectors of all the dislocation lines that cross the surface. It should 
also be noted here that the subsidiary constraint 



that the dislocation density be solenoidal must also be imposed to ensure that dislocation 
lines do not terminate within the solid, except perhaps at some internal interface such as 
a grain boundary. The necessity of incorporating this requirement was demonstrated by 
Nabarro. ^ 

Equation ([l|) is correct microscopically, but involves the unknown flux j^j , now a second 
rank tensor. We note that this equation has already been used as the starting point in other 
theories of dislocation dynamics [0,0] and that our approach differs in our choice of the flux 
jmj- As discussed earlier, expressions for the flux are derived invoking the Peach-Koehler 
force and continuum elasticity, and therefore yield equations that are invariant under time 
reversal. We suggest instead that dislocation motion is overdamped at the coarse-grained 
scale of motion, and that a constitutive equation for the flux can be introduced such that the 
dislocation density tensor evolves to minimize the local free energy of the defected material. 
The origin of this dissipative motion lies both in interactions with other microscopic degrees 
of freedom and interactions with other dislocations within the coarse-graining cell. As a 
consequence, our equations are manifestly not invariant under time reversal. 

As described earlier, the free energy of the crystal is a functional of the dislocation density 
tensor, and a non-singular strain ef^^*^, which describes phonon modes. It is useful to regard 




(2) 
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these two tensor fields as part of tlie basis of an abstract Hilbert space of slow variables 
which characterizes the system. Further, if these two fields are mutually orthogonal then, 
following Nelson and Toner H, 



NS 



NS 



+ Fd [Pij] , (4) 



where F/vg and Fd are the non-singular and the dislocation parts of the free energy, respec- 
tively. 

The spontaneous evolution of the system from some initial dislocation configuration 
towards equilibrium results in a decrease in its free energy (at constant temperature). We 
neglect in what follows the nonsingular contribution to free energy dissipation since for the 
phenomena of interest the associated time scale for phonon propagation is expected to be 
much shorter than that for dislocation motion. With this in mind, we have that the time 
rate of change of the free energy is given by 



dt J \^Pij J \ 9t 
where S/Spij stands for functional or variational derivative with respect to pij, the integral 
extends over the entire system, and we have chosen appropriate boundary conditions so that 
surface integrals vanish. Upon substituting the conservation law (Eq. (|l|)) and integrating by 
parts, the requirements that the total free energy be a monotonically decreasing function of 
time, and that there be a linear relation between forces and fluxes, leads to the identification 
of the dissipative dislocation density flux, 

• __R Al^^^ (f^^ 

where the fourth-rank tensor, B^jst plays the role of kinetic or Onsager coefficient, and 
embodies the complicated microscopic dynamics that allows the reduction in free energy at 
the coarse-grained scale. Substitution of Eq. into the conservation law, Eq. (p, leads 
to the dynamical equation for the dislocation density tensor. 
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Further simplification is possible when the kinetic attributes of the medium under con- 
sideration are spatially isotropic. In this case the Curie principle |]19| constrains the form of 
the tensor of kinetic coefficients to 



mjst 



B 



(7) 



where B is the single required kinetic coefficient. It is also possible to modify this de- 
scription in order to distinguish between glide and climb kinetics, though we will leave this 
modification for later work. As our formulation focuses on the coarse-grained dislocation 
density, the phenomenological constant B can only provide limited insight into the details 
of the mechanism for defect motion since, for example, the detailed geometry of slip planes 
and directions has been averaged over. Nevertheless, as will be shown below, one can ex- 
tract interesting microstructural information from our model, and so there is an implicit 
information trade-off in this approach to the problem. 

The only quantity that remains to be specified is the free energy density corresponding 
to the singular part of the deformation field, and to express it explicitly as a function of 
the dislocation density tensor. Kosevich has already calculated the elastic energy of 
interaction between sets of dislocations, and this result was later adopted by Nelson and 
Toner for equilibrium calculations [^. The result is most conveniently expressed in 
terms of an integration in reciprocal space by 



(8) 



with 



(9) 



where the main kernel K^j^i is given in terms of the shear modulus p and Poisson ration u 
by 



QikQjl + CilCkj + Y 



2iy 



— V 



■CijCki 



(10) 



and the core kernel Kf^j^i is given in terms of the edge and screw core energies (per unit 
length), Ee and Eg, by 

Ky{q) = 2EAkSji + 2{Es - E,)5,,5u. (11) 

It is advantageous to express these kernels in terms of the transverse and rotational g-space 
projection operators 

Q.. = - (12) 
Cij = (1^) 

respectively, as will be evident below. 

Before we proceed any further, we point out here that, upon substituting the densities 
corresponding to single isolated dislocations into Eq. (P), one recovers the familiar Peach- 
Koehler interaction 0. For example, consider two infinite, parallel screw dislocations, with 
Burgers vectors fo*-^-* and h^'^\ which are located at the origin and at a, respectively. The 
Fourier transforms of the dislocation densities corresponding to these defects are then 

p;j)(g-)=27r6«5,3 5,3%.), (14) 
pSf (g) = 27r 6(2)^^3 5^.3 5(g,)e-^- (15) 

and the line direction is along the 2;-axis. Neglecting the core contributions for the moment, 
one finds that the interaction free energy per unit dislocation length is 

implying that the associated force (per unit length) = fxb^^^b^'^^ /27ia, as expected. Thus, the 
main kernel given by Eq. (|T0]) embodies long-ranged interactions among dislocations. 

The effect of an external stress field with components a^j on a collection of dislocations 
can also be incorporated by the addition to the free energy F^, of a term representing the 
interaction of the dislocation density with the external field. Following Mura [Q, we write 
Fext as 
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F,,, = - j dhale:^, (17) 

where e*,- is the eigenstrain associated with a particular dislocation density tensor configu- 
ration, and the total free energy is then Ftot = + F^xt- We will elaborate further on this 
issue below. 

By combining the results summarized above, one finally arrives at the equation of motion 
for the dislocation density tensor, 

dt Spat{-q) 



III. STATIONARY SOLUTIONS 

Clearly, equation (0) is linear in the dislocation density tensor and will therefore only 
describe the linear relaxation of a defected system towards the asymptotic state of thermo- 
dynamic equilibrium, given by the minimization of Ffot- The relaxation rate as a function 
of wavevector q can be directly obtained from the eigenvalues of the kernel of Eq. (|1^). 

It is of interest to first study those eigenmodes of zero eigenvalue; i.e., the stationary 
solutions. We start by examining the equilibrium solutions, i.e., configurations of defects 
that minimize the free energy. 



A. Equilibrium configuration without external stress 

We first consider the free energy in the absence of external stresses and seek those 
densities that satisfy 



SF, 



kLM + KtpM] PpM = 0. (19) 



5pat{.-q) 

The kernel can be regarded as a 9x9 Voigt matrix with components Kap (1 < < 9) 
and, since Kap = KjSa, its eigenvalue spectrum is real. This kernel is non-singular and 
therefore the only eigenstates with zero eigenvalue are the trivial ones: ppq = 0, i.e., a 
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defect free medium. This fact inmediately follows from the core contribution to the free 
energy (Eq. (|11|)). However, Eq. (|T9]) also supports defected equilibrium states for a 
system constrained to having nonzero densities ppq. The Lagrange multipliers for such a 
constrained minimization can be chosen so as to cancel the core contribution, so that the 
problem therefore reduces to finding the (non-trivial) eigenstates of the kernel K^^^^ which 
have zero eigenvalue. This latter kernel is singular, hence one expects that there will be at 
least one such zero eigenstate. 

For this purpose it is convenient to express the dislocation density in terms of the op- 
erators which comprise the main kernel. For a rotationally invariant system one can form 
three second-rank tensors from q: 6ij, qiQj and tij^qk- We therefore use the decomposition 

m 



p,,{q) = 5,,X{q) + P,,Y{q) + Q,Z{q), (20) 

where the longitudinal projector Pij = Sij—Qij, and where X, Y and Z are three g-dependent 
amplitudes. Such a decomposition has the advantage that it permits the classification of 
the density in terms of its parity since, for example, Pij and Cij have different signatures 
under parity. Also, it is possible to determine the edge and screw character of the density 
by noting that for an edge (screw) dislocation the Burgers vector b is perpendicular (parallel 
or anti-parallel) to its line direction. Hence, pij is of edge (screw) character if pij ^ for 
^ 7^ J = j)- The effect of the main kernel on the dislocation density can then be determined 
from the results of the operator algebra summarized below: 

CsaCaq — Q sq (^-^) 
CsaQap Cgp QsaCap (^^) 
QiaQab — Qib (2^) 

With the preceding relations is mind, one can immediately identify several degenerate 
eigenstates with zero eigenvalue by observing that 
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0, 



(24) 



atpq^pg 



0, 



(25) 



and, consequently, Ppq is also an eigenstate. For the problem under consideration, then, 
Qpq is an acceptable eigenstate as it is solenoidal while 6pq clearly does not satisfy this 
requirement. In the former case, the transverse operator can be represented in real space by 
performing an inverse Fourier transformation to obtain 



Thus, the dislocation density falls off like and has a quadrupolar angular dependence. 

There are other degenerate eigenstates which correspond to more spatially localized 
dislocation densities and that break rotational invariance. As an example, consider the line 
along qi in reciprocal space 



which is an eigenstate of the main kernel with zero eigenvalue, given that A is a dimensional 
constant that is independent of wavevector. This particular combination of delta functions 
was chosen to ensure that the density is solenoidal. Now, this line in reciprocal space maps 
into the X2—X3 plane in real space. Further, the fact that the only two non-zero components 
are P22 and P33 implies that the density has screw character. In short, Eq. (^Tf) corresponds 
to a two-dimensional dislocation "wall" in real space. The screw character of the wall implies 
that it can be regarded, in some sense, as a prototypical (low-angle) twist grain boundary. 
It is interesting to note that such a regular distribution of dislocations might be expected 
from the tendency of dislocations to polygonize. In this model, however, such a structure 
has been obtained without postulating the existence of additional "chemical" forces, related 
to vacancy concentration at dislocations, which inhibit climb. Walls of other orientations 
are, of course, also possible. 




(26) 



Ppg = A [6pg - 6pl6gl] 5(^^)5(^5), 



(27) 



12 



B. With External Stress 



In many cases of interest the defected system will be subjected to an externally imposed 
stress field. Equation (|T^) also applies in this case, with the free energy explicitly including 
the contribution arising from the imposed stress. Now, in order to apply our previous results, 
it is first necessary to relate the eigenstrain in Eq. (|T7p to the dislocation density. Such 
a relation may be constructed by noting that the eigenstrains for this problem arise from 
plastic distortion, and that the essential connection between plastic strain and the Burgers 
vector of a dislocation can be expressed in terms of the closure failure of a Burgers circuit 
that encircles the dislocation. Generalizing this idea to a collection of dislocations leads to 

Piji^) = -(^ilm-Q^, (28) 

where the components of the plastic distortion tensor Pfj are related to the eigenstrain by 
e*j = (l/2)(/5^' + /3^). It should be noted here that Eq. (|28D by itself does not uniquely define 
the plastic strain. In particular, it does not specify the divergence of the plastic strain and 
so some additional input, such as the specification of a gauge, is required. The role of gauge 
invariance in the field theory of dislocations has been discussed elsewhere. [^. In the present 



context this additional input is contained implicitly in Eq. (^ since the dislocation flux is 
related to the time-rate change of the plastic distortion by jmk = —de^i^/dt, a relation 



indicated by Kroner and Rieder [Q. Further, from the form on the kinetic coefficient 
employed in Eq. (^), it is evident that jmm = 0, and so there is no time rate change of the 
volume associated with plastic deformation. 

With this in mind, we consider here, for the purpose of illustration, the specific case of 
a symmetric, solenoidal external stress tensor given in reciprocal space by 

aUq)=aiq)Qij, (29) 

where a{q) is a function of the wavenumber q. For this case F^^t becomes 

- / -^iQ)Q^^<A-^ = ^ / ^a,p.,(-g), (30) 
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ext 



where Eq. ( pSD has been used to relate the transverse projection of the dislocation density to 
the rotational projection. So, the state of equilibrium is defined by the eigenvalue equation 

= Kt,,i^PpM + -^c^t = 0. (31) 

'jpat[—q) q 

A solution to this eigenvalue problem can be obtained by noting that Cij itself is an eigenstate 



of the main kernel since 



KijkiCki = 2-^—Cij, (32) 
and so an equilibrium configuration (in reciprocal space) is given by 

or in real space by 

PP,{^ = J d'r'a{f')Lp,{f-r'), (34) 



where 



Lp^irl = Jil^'-IS^. (35) 



Thus, the application of a "transverse" stress field results in a dislocation density with 
"rotational" character. 

Since, as indicated above, Eq. (|]) relates the time-rate change of the eigenstrain to 
the dislocation density, it is advisable in the case of an arbitrary (symmetric) external 
stress to first develop an equation for the temporal evolution of the Fourier transform of the 
eigenstrain, e*j{q), and then determine the dislocation density by using the Fourier transform 
ofEq. (|2|) 

Pij(q) = -iqCiniP^j. (36) 

or the symmetrized form 

q"^ Gilm^jpk Plp^mk — Viji (37) 
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where rjij is Kroner's incompatibihty tensor p3[, which is essentially the Ricci tensor and 
therefore a measure of curvature. Such an equation can be determined by using Equa- 
tions d), (0), (|T9|) and to obtain 



de 



mk 



dt 



mkst 



QsqQtp ~l~ QspQtq 



2p 
1 - u 



QstQ 



PQ ' ^pq 



-pq 



(38) 



Thus, both the dislocation density and the eigenstrain can be determined from this for- 
mulation. Further, since the dependence of the eigenstrain on time is dependent on the 
dislocation flux, the value of the eigenstrain will depend upon the time history of the sys- 
tem, as expected. 

Equations ( ^6|) and (^) can be used to identify the dislocation density corresponding to 
a specific plastic strain. For example, in the case of no externally applied stress field, the 
plastic distortion Pfj = Pij corresponds to pij = 0, whereas the plastic distortion = Cij 
(and therefore efj =0) corresponds to pij = Qij. In each case, the plastic strain and the 
corresponding dislocation density of these equilibrium states are time-independent. 

Finally we address the issue of whether there exist other stationary solutions of the 
dynamic equations that are not minima of the free energy. We show next that the tensor 
CmiCsaBmjst, the generalized kinetic coefficient in Eq. (p!8|), commutes with the linear kernel 

for the class of solutions that do not break rotational invariance. Hence all the stationary 
solutions of the dynamics are minima of the free energy. Given the general decomposition for 
a rotationally invariant system of Eq. (|20D , and the fact that 6ij, Qij and Cij are eigenstates 
of K^, it is sufficient to show that 6ij, Qij and Cij are also eigenstates of CmiCsaBmjst- 

First, by using Eqs. (|^) and (|3), we find 



C miC saB mj stQ at 



BCr, 



2 i.^j'm ~l~ CfYij) ^^mjCg 



0. 



(39) 



The first term vanishes because Cij is an antisymmetric tensor. The second vanishes because 
Cij is proportional to the Levi-Civita tensor. Similarly, by using Eqs. ([7D, (pID and (|^) we 
find 



C miC saB mj stC at 



-BC„ 



1 



{Q 



mj 



Qjm) ^^mjQs 



B -a. 



:C 



V ] 1 



(40) 
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where we have used the fact that Qij is symmetric and that Q 



ss 



2. Finally, 




0. 



(41) 



IV. CONCLUSIONS 



We have considered a formal coarse-graining procedure over elements of volume that 
contain a large number of dislocation lines, and argued that microscopic interactions among 
the lines within each element of volume give rise to dissipative macroscopic dynamics. A 
constitutive law relating the dislocation density tensor and its associated flux has been 
introduced such that the evolution of the defected system is driven by the minimization of 
its plastic free energy. A standard procedure based on linear irreversible thermodynamics 
leads to the equations of motion for the dislocation density. 

For the case of an isotropic medium, the energy of interaction between dislocations in 
terms of the dislocation density tensor had already been calculated by Kosevich. We have 
further shown that in this case the tensor of kinetic coefficients that appears in the kinetic 
law requires only the specification of a scalar quantity, and that it commutes with the kernel 
of the linear operator driving energy minimization. Some stationary solutions have been 
explicitly constructed, both in systems with and without externally imposed stresses. 

We have restricted our attention in this paper to idealized systems wherein there is 
no mechanism for dislocation production, such as a Frank-Read source. [0] It is possible, 
however, to incorporate a stochastic "source" term in Eq. @ satisfying a fiuctuation- 
dissipation relation to model dislocation production. This would be useful in examining, 
for example, the change in mechanical properties that attend dislocation formation, such as 
work hardening. 

Finally, given that the approach discussed here is based on having partially averaged 
microscopic degrees of freedom, it seems of interest to relate it to results from microscopic 
descriptions such as the ones outlined in the Introduction. For example, one might consider 
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the dynamical evolution that results from these other models and then invoke a coarse- 
graining procedure. This, of course, would require systems that contain a large number of 
dislocations and, as mentioned above, algorithms for such large systems have or are currently 
being developed. Thus, it may be possible to bridge the disparate length scales between both 
descriptions. 
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